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Abstract 



A study of dijet production in proton-proton collisions was performed at = 7TeV 
for jets with pj > 35GeV and |y| < 4.7 using data collected with the CMS detec- 
tor at the LHC in 2010. Events with at least one pair of jets are denoted as "inclu- 
sive". Events with exactly one pair of jets are called "exclusive". The ratio of the 
cross section of all pairwise combinations of jets to the exclusive dijet cross section as 
a function of the rapidity difference between jets |Ay| is measured for the first time 
up to | Ay | = 9.2. The ratio of the cross section for the pair consisting of the most 
forward and the most backward jet from the inclusive sample to the exclusive dijet 
cross section is also presented. The predictions of the Monte Carlo event generators 
PYTHIA6 and PYTHIA8 agree with the measurements. In both ratios the HERWIG++ 
generator exhibits a more pronounced rise versus | Ay\ than observed in the data. The 
BFKL-motivated generators CASCADE and HEJ+ ARIADNE predict for these ratios a 
significantly stronger rise than observed. 
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The measurement of inclusive jet production in pp collisions provides an important testing 
ground for the Standard Model. Inclusive jet production is well described at LHC energies, 
over a wide range in jet transverse momentum and rapidity flU |2|, by calculations at next-to- 
leading-order (NLO) in perturbative quantum chromodynamics (QCD) using the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) approach ||3]-|7| and collinear factorization. The ra- 
pidity y is defined as y = (1/2) log[(£ + p z )/ (E — p z )\, where E is the jet energy and p z is 
the component of the jet momentum along the beam axis. As shown in [8], the production of 
dijets with invariant mass above 165 GeV and |y| less than 2.5 is also well described by NLO 
predictions. When jets are well separated in rapidity the description of the data becomes worse 
®- 

When the collision energy y/s is considerably larger than the hard scattering scale given by 
the jet transverse momentum, pj, the average number of produced jets grows rapidly, along 
with the phase space available in rapidity. This kinematic regime is expected to be described 
by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution [|9l4TT| and kj factorization fl2T[l4fl . 
An effective theory has been developed which describes strong interactions in this kinematic 
domain IITBll . 

The ratio of the dijet production cross section in "inclusive" events to that in "exclusive" events, 
R mcl = cr mc1 / (7 excl , as a function of the rapidity separation | Ay| between two jets, is a sensitive 
probe of effects beyond collinear factorization [TBI . Only jets with transverse momenta above 
a minimal value of p™ m are considered. Events with at least one pair of jets are denoted as 
"inclusive". Events with exactly one pair of jets are called "exclusive". In the inclusive case, 
the rapidity separation is evaluated for each pairwise combination of jets above threshold [p~6|. 
Mueller-Navelet jet pairs iTPTll are a subset of the inclusive dijet class. In this case only the jet at 
highest rapidity (i.e. most forward) and that at lowest rapidity (most backward) are considered. 
At low \Ay\, the inclusive dijet ratio, R , is always larger than the corresponding Mueller- 
Navelet dijet ratio, R MN = c MN / cj excl , because all the jets in the rapidity interval between the 
Mueller-Navelet jets contribute to the inclusive cross section. 

From the theoretical point of view, an advantage of the ratios R mcl and R MN with respect to the 
individual dijet production cross sections is that the influence of the uncertainty of the parton 
distribution functions is greatly reduced; in addition, the ratios are particularly sensitive to 
the parton radiation pattern [16]. At large enough energies, the parton subprocesses involve a 
large number of partons with comparable transverse energies. Such subprocesses, governed by 
BFKL evolution, lead to an increase of the ratios with increasing \Ay\ ffl6T[T8| . Experimentally, 
an additional advantage of the R mcl and R MN ratios defined above is the cancellation of most 
of the systematic uncertainties affecting both numerator and denominator. 

Earlier measurements of Mueller-Navelet jets with \Ay\ up to 5 were made at the Tevatron by 
the DO experiment ffl9l[20 ]. DO did not find indications of BFKL effects in the azimuthal decor- 
relation data [19J; however, a stronger than expected \/s dependence of the dijet production 
cross section was observed H20| when comparing the data at \Js = 630 and 1800 GeV for dijets 
with large rapidity separation. The ATLAS Collaboration recently studied various dijet produc- 
tion cross section ratios at \/s = 7TeV for \Ay\ < 6 |2T|; the results are roughly in agreement 
with the DGLAP predictions. At the HERA ep collider, deviations from DGLAP were observed 
in studies of forward jets by the ZEUS |22l |23| and HI [24, 25 [ Collaborations. However, no 
compelling evidence for BFKL effects was found. 

The component of the CMS detector [26 1 most relevant for this analysis is the calorimeter sys- 
tem extending to pseudorapidities \rj\ = 5.2, where rj = — log[tan(#/2)], and 9 is the polar 
angle relative to the anticlockwise proton beam direction. The crystal electromagnetic calori- 
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meter (ECAL) and the brass/ scintillator hadronic calorimeter (HCAL) extend to pseudorapidi- 
ties \rj | = 3.0. The HCAL cells map to an array of ECAL crystals to form calorimeter towers 
projecting radially outwards from the nominal interaction point. The pseudorapidity region 
3.0 < \rj I < 5.2 is covered by the hadronic forward (HF) calorimeter, which consists of steel ab- 
sorber wedges with embedded radiation-hard quartz fibers, oriented parallel to the beam direc- 
tion. The calorimeter towers in the barrel region have segmentation of Arj x A(p = 0.087 x 0.087, 
becoming progressively larger in the endcap and forward regions (A// x Acp = 0.175 x 0.175 at 
rj ~ 4.5). 

The CMS trigger system consists of a hardware Level-1 trigger and a software high-level trig- 
ger. Jets formed online by the trigger system use ECAL, HCAL and HF inputs for energy 
clustering and are not corrected for the jet energy response. The minimum-bias trigger was 
defined as the coincidence of a signal from one of the two beam scintillation counters covering 
the range 3.23 < \rj\ < 4.65 and a signal from one of the beam pick-up timing devices. 

The data were collected in 2010, when LHC collided protons at s/s = 7TeV. Various triggers 
were used to cover the largest possible range in rapidity separation between jets. Dijets with 
moderate | Ay\ were selected by a single-jet trigger with a threshold on the raw (uncorrected) jet 
transverse momentum of 15 GeV. This trigger was significantly prescaled as the instantaneous 
luminosity increased, and the effective integrated luminosity recorded with it is ~ 33nb _1 . 

Dijet events with large rapidity separation are rare. Therefore, a dedicated trigger for forward- 
backward dijets was developed. This forward-backward-dijet trigger selects events with two 
jets in opposite hemispheres and \f]\ > 3.0, and jet raw transverse momentum pj > 15 GeV. It 
was operated with moderate prescaling, and the effective integrated luminosity recorded with 
it is ~ 5pb _1 . This allowed to collect a number of dijet events at high \Ay\ values more than 
100 times larger than with the single-jet trigger. 

The trigger efficiency was measured by means of a control sample selected with the minimum- 
bias trigger. The single-jet trigger was found to be 100% efficient for dijets with corrected 
px > 35 GeV. The single-jet trigger was also used for the determination of the efficiency of 
the forward-backward-dijet trigger. The latter was 100% efficient for dijets with pj > 35 GeV. 

Jets were reconstructed offline from the energy depositions in the calorimeter towers, clustered 
with the anti-kj algorithm |27l [28] with a distance parameter R = 0.5. In the reconstruction 
process, the contribution from each tower was assigned a momentum, the absolute value and 
the direction of which were given by the energy measured in the tower, and the coordinates of 
the tower, respectively. The raw jet energy was obtained from the sum of the tower energies, 
and the raw jet momentum from the vectorial sum of the tower momenta. The raw jet ener- 
gies were then corrected to establish a uniform relative response of the calorimeter in rj and 
a calibrated absolute response in transverse momentum pj II29II . The jet energy resolution for 
calorimeter jets with p T ~ 35 GeV is about 22% for \t]\ < 0.5 and about 10% for 4 < \t]\ < 4.5 
l[3"0") . The uncertainty of the jet energy calibration for jets with pj ~ 35 GeV depends on rj and 
is ~ 7 - 8% EU. 

In the offline analysis, at least one well-reconstructed primary vertex is required to be present 
within ±24 cm of the nominal interaction point along the beamline |3"T|. In order to reduce the 
sensitivity to overlapping pp collisions (the so-called "pile-up" events), events with only one 
primary vertex reconstructed within the luminous region were used for the measurement. 

Loose jet quality cuts [32] were used to suppress the effect of calorimeter noise. Events with at 
least two jets with pj > 35 GeV and \y\ < 4.7 were selected; only jets satisfying these criteria 
were used for the analysis. All pairwise combinations of jets from the selected events entered 
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the inclusive distribution. For studies of Mueller-Navelet jets, only the pair consisting of the 
most forward and the most backward jet was considered. The exclusive dijet sample is a subset 
of the inclusive and Mueller-Navelet samples, and consists of events where exactly one pair 
of jets is found. The measured observables, R mcl and R MN , are defined as the ratios of the 
yield of inclusive or Mueller-Navelet dijets to the yield of exclusive dijets in a specific | Ay| bin, 
respectively. 

Detector effects were accounted for by applying bin-by-bin corrections derived from Monte 
Carlo (MC) simulations. Simulated events produced with the generators PYTHIA6 (version 
6.422) tune Z2 [33H35H and HERWIG++ ESI (version 2.4.2 with default settings) were passed 
through the full CMS detector simulation based on the Geant4 package [37], and were input 
to the same event-reconstruction program as used for the data. To quantify detector effects, 
the distributions obtained from detector-level quantities were compared to the distributions 
obtained at the level of stable particles (lifetime T such that ct > 10 mm). The ratios of the 
stable-particle level and detector-level quantities in a given bin were used to correct the data. 

Because of the finite px and y resolutions, the detector-level distributions deviate from the 
corresponding ones at stable-particle level. Events can migrate to and from the exclusive or the 
inclusive samples because of fluctuations of the measured transverse jet momentum around the 
p!p n threshold. The amount of these migrations was estimated not to exceed 20%. Similarly, the 
| Ay | value measured at detector level may fall into a different \Ay\ bin compared to the stable- 
particle level; this effect, for the present data, is typically around 5-10% and reaches 15-25% at 
most. The influence of these migrations on the measured ratios R incl and R MN is minimal as 
the effects for numerator and denominator are similar. The uncorrected R mcl and R MN ratios 
are reasonably well reproduced by the PYTHIA6 events and less well by HERWIG++ events, 
both of which have been passed through the detector simulation. The px arid rj distributions of 
the jets at the detector level for both MC generators agree with the data within the jet energy 
scale uncertainty (see also [2]). The correction factors were therefore obtained with PYTHIA6, 
while their model dependence was estimated from the difference between the PYTHIA6 and 
HERWIG++ corrections. 

The following sources of systematic effects were considered: 

1. Uncertainty of the jet energy calibration. The uncertainty of the measurement was es- 
timated by shifting the jet energy scale (JES) by the pj- and ^-dependent uncertainties 
derived in 1291 . The resulting variation of the measurements does not exceed 4.2% for 
R M and 3.8% for R MN . 

2. Uncertainty of the corrections for detector effects. 

(a) Uncertainty due to model dependence of the correction factors. As discussed above, 
correction factors were determined by using HERWIG++ and PYTHIA6, and corrected 
measurements were obtained for each case. The difference in the results was taken 
as a measure of the model dependence of the correction factors; it does not exceed 
3.4% for R M and 3.3% for R MN . 

(b) Uncertainty related to the quality of the MC description of the jet px and y resolu- 
tions. These resolutions were modified by ±10% in the simulation as recommended 
in [|30|. A conservative estimate of the uncertainty associated to this effect does not 
exceed ±1.0%. 

3. An additional systematic effect may originate from the extra energy and jets due to pile- 
up collisions. As noted earlier, the measurement was restricted to events with only one 
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reconstructed primary vertex to reduce the impact of pile-up. However, a contribution 
from interactions undetected because of vertex reconstruction inefficiency may still be 
present. By comparing data taken at different instantaneous luminosities, an upper limit 
of 1.3% was estimated for this effect. The value was obtained for pairs of jets in the 
forward region, where the impact of pile-up on jet reconstruction is more significant [29 1; 
the effect is thus expected to be smaller for jets at central rapidities. 

The total systematic uncertainty was calculated as the quadratic sum of the individual uncer- 
tainties. The total uncertainty is \Ay\ dependent but always smaller than 5.6% for R mcl and 
4.8% for R MN (Table 1). 

Table 1: Sources of systematic effects and associated uncertainties. The ranges correspond to 
the variation of the uncertainty with \Ay\. For different uncertainty sources, the minimum 
and maximum values may correspond to different |Ay| bins. For asymmetric uncertainties the 
upper and lower limits are shown. 



source 


R inci uncertainty (%) 


r mn uncertainty (%) 


jet energy scale 
uncertainty of detector corrections 
pile-up 


+ (2.2-4.2) 
-(1.0-3.0) 

±(1.5-3.5) 
< 1.3 


+ (0.2-3.8) 
-(0.2-2.3) 

±(0.1-3.4) 
< 1.3 


total 


+ (2.8-5.6) 
-(2.7-4.5) 


+ (0.2-4.8) 
-(0.2-3.7) 



The measured ratios, corrected for detector effects, were compared to the predictions of sev- 
eral MC generators at the stable-particle level. A DGLAP leading-order parton shower ap- 
proach is used in PYTHIA6 (version 6.422) tune Z2 (35l, PYTHIA8 (version 8.145) [38] tune 
4C EU and HERWIG++ (version 2.5.1) tune UE-7000-EE-3 ESUQl For the simulation of the 
non-perturbative fragmentation, PYTHIA uses the Lund string model, and HERWIG++ the clus- 
ter fragmentation approach. The tunes mentioned above include multiple parton interactions 
(MPI) as a part of the underlying event (UE) modeling. The dijet observables might be affected 
by MPI through jets that do not originate from the same hard interaction. The effect was esti- 
mated by switching off the MPI options in PYTHIA6 and HERWIG++; no significant change in 
the predictions was observed. 

The Monte Carlo generators CASCADE (version 2.2.03) EH and hej (version 1.3.2) 021 are mo- 
tivated by the leading-logarithmic BFKL approach and incorporate parts of a next-to-leading 
logarithmic approximation. The HEJ generator produces parton-level jets; the corresponding 
showers were produced with the ARIADNE program [43 1. The HEJ+ ARIADNE package [44j ver- 
sion 0.99b, consisting of hej 1.3.2 and ARIADNE 4.12, was used. 

The ratio R mcl of inclusive to exclusive dijet production as a function of \Ay\ is presented in 
Fig. [I] (left panel). On average the inclusive cross section is 1.2 - 1.5 times larger than the 
exclusive cross section. The ratio R mcl grows with increasing \Ay\, as expected because of the 
larger phase space for hard parton radiation. At the highest \Ay\, R mcl is expected to decrease 
because energy-momentum conservation suppresses the emission of extra jets. The \Ay\ value 
where R mcl starts to decrease varies from one MC generator to another, as can be seen in Fig. [l] 

The predictions from PYTHIA6 and PYTHIA8 agree with the measurement. HERWIG++ overes- 
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Figure 1: Ratios of the inclusive to exclusive dijet cross sections as a function of the rapidity 
separation \Ay\ between the two jets, R mcl (left panel) and R MN (right panel), compared to the 
predictions of the DGLAP-based MC generators PYTHIA6, PYTHIA8 and HERWIG++, as well 
as of CASCADE and HEJ+ ARIADNE, which incorporate elements of the BFKL approach. The 
shaded band indicates the size of the total systematic uncertainty of the data. Statistical uncer- 
tainties are smaller than the symbol sizes. Because of limitations in the CASCADE generator, it 
was not possible to obtain a reliable prediction for |Ay| > 8. 



co 1-5 

t& 

^1.4 
O 

2 1.3 
1.2 
1.1 



CMS, pp,\fs = 7TeV 



0.8. 







: ■ data 


dijets, R incl : 


\ PYTHIA6Z2 


P T > 35 GeV - 


- PYTHIA8 4C 


|y|<4.7 I 


; HERWIG++ 














" 





1 



8 9 

|Ay| 



co 1-5 
to 

^1.4 
O 

1.2 
1.1 
1 

0.9 
0.8. 



data 

PYTHIA6 Z2 
PYTHIA8 4C 
HERWIG++ 



CMS, pp,Vs = 7TeV 

I i i i i I i i 

dijets, R MN 
p T > 35 GeV - 
|y|<4.7 



1 



8 9 

|Ay| 



Figure 2: Predictions for R incl (left) and R MN (right) from DGLAP-based MC generators pre- 
sented as ratio to data corrected for detector effects. Both BFKL-motivated generators CASCADE 
and HEJ+ ARIADNE (not shown) lead to a MC/ data ratio well above unity. The shaded band 
indicates the size of the total systematic uncertainty of the data while statistical uncertainties 
are shown as bars. 
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timates the ratio R mcl at medium and large rapidity intervals. A detailed comparison between 
the data and the predictions of the DGLAP-based MC generators is presented as a ratio in Fig. [2] 
(left panel). It was checked explicitly that the results obtained from PYTHIA6 and HERWIG++ at 
parton level are close to the corresponding ones at stable-particle level. The different behaviour 
of PYTHIA6 and HERWIG++ is also observed at parton level. 

The ratio R MN and the corresponding MC to data ratio are presented in the right panels of Figs. 
[l]-[2j At large |Ay|, R MN approaches R mcl as extra jet radiation contributing to R mcl tends to con- 
centrate at moderate rapidities. The quality of the predictions of the DGLAP-based MC gen- 
erators for R^™ is similar to those for R incl . The MC generators CASCADE and HEJ+ ARIADNE 
considerably overestimate the measurements of both R mcl and R MN . 

The ATLAS Collaboration measured several observables for dijets as functions of the rapidity 
separation between jets |2l]. One is the fraction of events without jets with pj above pj et0 in 
the rapidity interval between the most forward and the most backward jets, which is equal 
to 1/R MN . A difference between data and the parton-level HEJ prediction at large |Ay| for 
70 < (p T ) < 90GeV, p\ eto = (p T ) and |Ay| < 6 was reported (Fig. 8 from 121]), which is in 
qualitative agreement with the result presented here. 

To conclude, the first measurement of the ratios R mcl and R MN in a wide range of rapidity sep- 
aration, up to | Ay | = 9.2, in proton-proton collisions at >/s = 7TeV was presented. A moderate 
rise of the ratio of the inclusive to exclusive dijet production cross sections as a function of | Ay | 
is observed. The predictions of the PYTHIA6 and PYTHIA8 generators agree with the measure- 
ments. The predictions of the HERWIG++ generator are larger than the measurement especially 
at large |Ay|. The BFKL-motivated generators CASCADE and HEJ+ ARIADNE predict for these 
ratios a significantly stronger rise than observed. The moderate rise of the measured dijet ra- 
tios indicates that the BFKL effects are not dominant for jets with pj > 35 GeV at the present 
collision energy of 7 TeV. 
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